Discovery of very high energy gamma-rays from the flat 
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Abstract. 3C 279 is one of the best studied flat spectrum radio quasars located at a comparatively large redshift of z = 0.536. 
Observations in the very high energy band of such distant sources were impossible until recently due to the expected steep 
energy spectrum and the strong gamma-ray attenuation by the extragalactic background light photon field, which conspire 
to make the source visible only with a low energy threshold. Here the detection of a significant gamma-ray signal from 
3C 279 at very high energies (E > 75 GeV) during a flare in early 2006 is reported. Implications of its energy spectrum on the 
current understanding of the extragalactic background light and very high energy gamma-ray emission mechanism models are 
discussed. 
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INTRODUCTION 

The flat spectrum radio quasar (FSRQ) 3C 279, located 
at a redshift of z = 0.536, is one of the brightest extra- 
galactic objects in the gamma-ray sky 1^1 and it is an ex- 
ceptionally variable source at various energy bands, in- 
cluding very strong gamma-ray flares recorded by satel- 
lite detector EGRET in 1991 and 1996 U. Its detection 
at very high energy (VHE, defined here as £ > 50 GeV) 
gamma-rays by the MAGIC telescope "T] is the first 
detection of an FSRQ in this energy band as well as the 
farthest object detected at these wavelengths, as the high- 
est redshift previously observed was z — 0.212 |5]. 

The EGRET experiment detected 66 blazars: 5 1 FSRQ 
and 15 BL Lac objects [6]. Blazars are thought to be 
powered by accretion of matter onto supermassive black 
holes residing in the centers of galaxies and ejecting rel- 
ativistic jets with small angles to the line of sight 0]. 
Except for the FSRQ 3C 279, all the blazars detected at 
VHE up to know are of the BL Lac type, being the dif- 
ference between this two subclasses that FSRQs show 
optical emission lines while in BL Lacs these are very 
weak or absent. The spectral energy distribution (SED) 
of blazars is characterized by two broad bumps, the first 
one peaking between the infrared and the X-ray band and 
the second one in the gamma-ray domain. The first peak 
is usually attributed to synchrotron emission from rela- 



tivistic electrons, while the high-energy peak is presum- 
ably due to the inverse Compton scattering of low-energy 
photons by these electrons. The soft photon population 
can be the synchrotron photons themselves (synchrotron 
self-Compton or SSC mechanism f§\) or ambient radia- 
tion (external Compton or EC mechanism [9, 10]). Other 
scenarios (hadronic models) involve gamma-ray emis- 
sion due to accelerated protons and ions iflTlll^ . 

Apart from being the first FSRQ detected in the VHE 
band, the detection of 3C 279 is important because 
gamma-rays from distant sources are expected to be 
strongly attenuated in intergalactic space by the possible 
interaction with low-energy photons (7+ 7 ^ e+ + e^). 
These photons form the extragalactic background light 
(EBL) lfl3ll and have been radiated by stars and galaxies 
in the course of cosmic history. The attenuation by the 
EBL increases with increasing distance to the source and 
energy of the emitted gamma-ray, therefore, VHE radia- 
tion from distant sources is strongly suppressed at high 
energies (for E > 300 GeV and z = 0.536 more than 90% 
of the emitted gamma-rays are expected to be absorbed). 

OBSERVATIONS AND ANALYSIS 

Observations of 3C 279 were performed using the 
MAGIC telescope (see ifTill for a detailed description). 



which is located on the Canary Island of La Palma 
and has a 17m-diameter tessellated reflector dish and 
a photomultiplier camera covering 3.5° field of view. It 
has a trigger threshold of 60 Ge V, an angular resolution 
of approximately 0.1° and an energy resolution above 
150 GeV of about 25%. 

Data were taken in on/off mode, which includes on 
events taken with the telescope pointing at the source 
and a similar amount of off events taken at a sky loca- 
tion near the source with very similar operating condi- 
tions. 3C 279 was observed over 10 nights between Jan- 
uary and April 2006 for a total of 9.7 hours. Simulta- 
neous optical R-band observations with the 1 .03-m tele- 
scope at the Tuorla Observatory (Finland) and the 35-cm 
KVA telescope on La Palma revealed that during MAGIC 
observations 3C 279 was in high optical state, a factor 
of 2 above its long-term baseline flux (host galaxy sub- 
tracted). 

The data anaysis was carried out using the standard 
MAGIC analysis and reconstruction software |14]. To 
search for a gamma-ray signal, first the standard calibra- 
tion of the data was performed. Then an image clean- 
ing procedure was applied using the amplitude of the 
calibrated signals. The next step included the parame- 
terization of the shower images |15]. Hadronic back- 
ground suppression was achieved using the Random For- 
est (RF) method lfT6l ITtIi . where for each event the so- 
called hadronness (/i) parameter is computed based on 
the image parameters. Moreover, the RF method was also 
used for the energy estimation using a Monte Carlo sim- 
ulated gamma-ray sample with the same zenith angle dis- 
tribution as the data sample. All image parameters were 
checked for consistency between on and off data. 



RESULTS 

The excess events were obtained by subtracting suitably 
normalized off from the on data. Cuts in h and alpha 
(the latter describes the direction of the main axis of the 
image) were optimized using Crab Nebula data taken at 
similar zenith angles. For the detection and light curve 
the cuts /; < 0.12 and alpha < 12° where used, resulting 
in a combined efficiency of 40% for gammas. For these 
cuts and a quality cut that removed small images the 
resulting energy threshold is 170 GeV. 

3C 279 was clearly detected on the night of 23 Febru- 
ary 2006 (Fig. [T]l with a significance of 6.2(7, which 
translates into 5.8(7 when corrected with a trial factor 
corresponding to 10 nights of observation. A less signif- 
icant excess on the previous night is visible in the light 
curve (Fig. 121). As determined by the test, the probalil- 
ity that the gamma-ray flux on all 10 nights was zero is 
2.3 X 10^^, corresponding to 5.0a in a Gaussian distri- 
bution. 
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FIGURE 1. Histogram of the image parameter alpha ob- 
tained for 3C 279 in the night of 23 February 2006 for on 
(dots) and off (crosses and filled area) pointings. The dotted 
line is a simple parabolic fit without linear term to the off events 
between 0° and 80°, serving to smoothen the background distri- 
bution for a better extrapolation towards alpha = 0. The num- 
ber of excess events is obtained with respect to this line. 
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FIGURE 2. MAGIC (top) and optical R-band (bottom) light 
curves obtained for 3C 279 from February to March 2006. The 
long-term baseline for the optical flux is 3 mJy (host galaxy 
subtracted). 



The observed differential VHE spectrum (Fig. [3]) was 
obtained with looser cuts to increase the number of 
gamma-ray candidates, and can be described by a power 
law with spectral index of O! = 4. 1 ± 0.7stat i 0.2syst. The 
measured integrated flux above 100 GeV on 23 February 



is 5.15±0.82stat±1.50, 



syst, 



X lO^^'cm-^s-i. 



DISCUSSION 

The EBL influences the propagation of gamma-ray emit- 
ted from distant sources resulting in an exponential de- 
crease with energy and a cutoff in the measured gamma- 
ray spectrum. Several models have been proposed for 
the EBL L 1 3i1 . All have limited predictive power for the 
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FIGURE 3. Differential energy spectrum of 3C 279 mea- 
sured by MAGIC. The grey area includes the combined sta- 
tistical (1(7) and systematic errors, and underlines the marginal 
significance of detections at high energy. The dotted line shows 
compatibility of the measured spectrum with a power law of 
photon index a = 4.1. The triangles are measurements cor- 
rected on the basis of the two models for EBL density discussed 
in the text. 
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FIGURE 4. EBL models at z = and direct measurements 
of the EBL photon density at various wavelengths. The upper 
limit derived from the detection of 3C 279 (black line) is based 
on 1 22] with parameters adapted to the results from the latest 
galaxy counts (see Supporting Material in 0] for details). The 
dotted line is an upper limit from |2?]. The shaded region 
corresponds to the range of frequencies where the MAGIC 
measured spectrum is sensitive. 

EBL density, particularly as a function of the redshift, as 
many details about star and galaxy evolution remain un- 
certain. To illustrate this uncertainty two extreme mod- 
els are used to calculate the source intrinsic spectrum: a 
model by Primack et al. LI 8.1 close to the lower limits set 
by galaxy counts |[l 9ll20ll and the 'fast-evolution' model 
by Stecker et al. [21] corresponding to the highest possi- 
ble attenuation. These models will be referred as low and 
high respectively and are shown in Figure |4l 

A power-law fit to the EBL-corrected spectral points 
results in an intrinsic photon index a* — 2.9 ± 0.9stat ± 
O.Ssyst (low absorption) and a* = 0.5 ± 1.2stat ±0.5syst 
(high absorption). The high model leads to an intrinsic 
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FIGURE 5. The gamma-ray horizon. The redshift region 
over which the gamma-ray horizon can be constrained by ob- 
servations has been extended up to z = 0.536. The prediction 
range of EBL models is illustrated by |18] (thick solid line) 
and [21] (dashed-dotted line). The tuned model of [ 22] (dashed 
line) represents an upper EBL limit based on our 3C 279 data, 
obtained on the assumption that the intrinsic photon index is 
> 1.5 (arrow corresponding to 3C 279). Limits obtained for 
other sources are shown by black arrows, most of which lie 
very close to the model |22]. The narrow filled band is the 
region allowed between this model and a maximum possible 
transparency (i.e., minimum EBL level) given by 1 18], which is 
nearly coincident with galaxy counts. The gray area indicates 
an optical depth T > 1, where the flux of gamma rays is strongly 
suppressed. To illustrate the strength of the attenuation in this 
area energies for T = 2 and T = 5 (thin black lines) are shown, 
again with 81811 as model. 



VHE spectrum difficult to reconcile with an extrapola- 
tion of the EGRET data and with general constraints in 
the spectral energy distribution. On the other hand, the 
low model apparently gives an acceptable result. The dis- 
tance at which the flux of photons of a given energy is 
attenuated by a factor e (the path corresponding to an op- 
tical depth T = 1) is called the gamma-ray horizon f24'l 
and is usually expressed as a function of the redshift. This 
energy-redshift relation is shown in Figure|5] 

Assuming a hardest intrinsic pho ton index a* = 1 .5 a 
model based on Kneiske et al. 1q3l Izsll can be tuned to 
give an upper limit to the EBL density in the 0.2 — 2 jtim 
range as shown in Fig.|4] This result supports the conclu- 
sion drawn from VHE measurements at lower redshifts 
^2^ that the EBL density is close to the lower limits set 
by galaxy counts. It is important to point out that this 
limit is only valid under the assumption that the intrin- 
sic photon index cannot be harder than 1.5, which is the 
hardest value given for EGRET sources and is also the 
hardest that can be obtained with classical leptonic emis- 
sion mechanisms. Alternative scenarios can produce in- 
trinsic spectra with a* < 1.5. SSC models with a narrow 
electron distribution can produce spectra with a* ~ 0.7 
ll27ll . and internal absorption by soft photons can also 
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FIGURE 6. Spectral energy distribution for 3C 279 at differ- 
ent source intensities, over tfie years from 1991 to 2003, taken 
from 1 34, 35] with MAGIC points (2006) at the high-energy 
end. The MAGIC points are corrected according to the two 
extreme EBL models discussed on the text. The spectrum cor- 
rected with the low EBL model gives clearly a more plausible 
result. An emission model and its individual components are 
also shown: the solid line represents the total model emission, 
fitted only to the MAGIC points unfolded with a low EBL den- 
sity and the simultaneous optical point from the KVA telescope. 
The individual components of the emission model are syn- 
chrotron radiation (dotted line), disk emission (long-dashed), 
synchrotron-self-Compton (short-dashed) and external Comp- 
ton (dot-dashed). 

harden the spectrum ll28i l29i [30l [sill . However, an ac- 
curate modehng of the internal absorption in 3C 279 
finds no important hardening of its spectrum in the en- 
ergy band measured by MAGIC HI]. 

The emission mechanism responsible for the observed 
VHE radiation remains uncertain. The spectrum cor- 
rected with a low level EBL can be reproduced by a 
one-zone SSCh-EC model |33] (Fig. |6]l while the spec- 
trum corrected with a high absorption can not be fitted 
by generic leptonic emission models. Studies including 
quasi-simultaneous optical and X-ray data |36] conclude 
that neither one-zone SSC nor one-zone EC models can 
reproduce the observed SED iIstIi . 
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